Search for: All records

We explore the relation between stellar surface density and gas surface density (the star–gas, or S-G, correlation) in a 20,000M_⊙simulation from the STAR FORmation in Gaseous Environments (starforge) project. We create synthetic observations based on the Spitzer and Herschel telescopes by modeling contamination by active galactic nuclei, smoothing based on angular resolution, cropping the field of view, and removing close neighbors and low-mass sources. We extract S-G properties such as the dense gas-mass fraction, the Class II:I ratio, and the S-G correlation (Σ_YSO/Σ_gas) from the simulation and compare them to observations of giant molecular clouds, young clusters, and star-forming regions, as well as to analytical models. We find that the simulation reproduces trends in the counts of young stellar objects and the median slope of the S-G correlation. This implies that the S-G correlation is not simply the result of observational biases, but is in fact a real effect. However, other statistics, such as the Class II:I ratio and dense gas-mass fraction, do not always match observed equivalents in nearby clouds. This motivates further observations covering the full simulation age range and more realistic modeling of cloud formation.

 

Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the starforge simulation suite – a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use nbody6++gpu to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyse the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e. global cloud surface densities below 0.15 g cm−2, star formation efficiencies below 15 per cent, and gas expulsion time-scales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have ∼1000 bound members, form in the densest and lowest velocity dispersion clouds, representing ∼32 and 39 per cent of the stars in the simulations, respectively. The cloud’s early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power-law relation at 10 Myr with a power-law index of 0.21. This correlation result in a distinct mass–size relationship for bound clusters. We also provide valuable constraints on the gas dispersal time-scale during the star formation process and analyse the implications for the formation of bound systems.

 

Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud, the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation underpredicts these values. The discrepancy is likely due to the lack of disc fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $\gt 1\, \mathrm{M}_{\rm \odot }$ stars. However, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anticorrelates with stellar density. This relationship also explains most multiplicity variations between runs, i.e. variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) also act to decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties.

 

One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies and the Central Molecular Zone. In this work, we present simulations from the STARFORGE project that are the first cloud-scale radiation-magnetohydrodynamic simulations that follow individual stars and include all relevant physical processes. The simulations include detailed gas thermodynamics, as well as stellar feedback in the form of protostellar jets, stellar radiation, winds, and supernovae. In this work, we focus on how stellar radiation, winds, and supernovae impact star-forming clouds. Radiative feedback plays a major role in quenching star formation and disrupting the cloud; however, the IMF peak is predominantly set by protostellar jet physics. We find that the effect of stellar winds is minor, and supernovae ‘occur too late’ to affect the IMF or quench star formation. We also investigate the effects of initial conditions on the IMF. We find that the IMF is insensitive to the initial turbulence, cloud mass, and cloud surface density, even though these parameters significantly shape the star formation history of the cloud, including the final star formation efficiency. Meanwhile, the characteristic stellar mass depends weakly on metallicity and the interstellar radiation field, which essentially set the average gas temperature. Finally, while turbulent driving and the level of magnetization strongly influence the star formation history, they only influence the high-mass slope of the IMF.

 

Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of individual, low-mass stars. We find that individual star-formation sites merge to form ever larger structures, while still accreting gas. Thus clusters are assembled through a series of mergers. During the cluster assembly process, a small fraction of stars are ejected from their clusters; we find no significant difference between the mass distribution of the ejected stellar population and that of stars inside clusters. The star-formation sites that are the building blocks of clusters start out mass segregated with one or a few massive stars at their centre. As they merge the newly formed clusters maintain this feature, causing them to have mass-segregated substructures without themselves being centrally condensed. The merged clusters relax to a centrally condensed mass-segregated configuration through dynamical interactions between their members, but this process does not finish before feedback expels the remaining gas from the cluster. In the simulated runs, the gas-free clusters then become unbound and breakup. We find that turbulent driving and a periodic cloud geometry can significantly reduce clustering and prevent gas expulsion. Meanwhile, the initial surface density and level of turbulence have little qualitative effect on cluster evolution, despite the significantly different star formation histories.

 

Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or ‘stirred’ turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called turbsphere) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a ∼k−2 velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results similar to previous simple spherical cloud setups, but external driving can suppress star formation considerably in the new setup. Periodic box simulations with the same dimensions and turbulence parameters form stars significantly slower, highlighting the importance of boundary conditions and the presence or absence of a global collapse mode in the results of star formation calculations.

 

We analyse the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $\sim 8{{\ \rm per\ cent}}$ star formation efficiency (SFE), and the first supernova at $8.3\, \rm Myr$ comes too late to influence star formation significantly. The per-free-fall SFE is dynamic, accelerating from 0 per cent to $\sim 18{{\ \rm per\ cent}}$ before dropping quickly to <1 per cent, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($\sim 1\, \mathrm{Myr}$) compact ($\sim 1\, \rm pc$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope α = −2 and a maximum mass of 55 M⊙. Stellar accretion takes $\sim 400\, \rm kyr$ on average, but $\gtrsim 1\,\rm Myr$ for >10 M⊙ stars, so massive stars finish growing latest. The fraction of stars in multiples increase as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties.

 

We characterize protostellar multiplicity in²⁰

Current address: Niels Bohr Institute, University of Copenhagen, Øster Voldgade 5–7, DK-1350, Copenhagen K, Denmark.

the Orion molecular clouds using Atacama Large Millimeter/submillimeter Array 0.87 mm and Very Large Array 9 mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as ∼20 au, and we consider source separations up to 10⁴au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30 ± 0.03 and 0.44 ± 0.03, respectively, considering separations from 20 to 10⁴au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at ∼100 au and ∼10³au, we argue that multiples with separations <500 au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at ≳10³au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10³au in regions of high young stellar object density. This finding may be evidence for the migration of companions from >10³au to <10³au, and that some companions between 10³and 10⁴au must be (or become) unbound.